Robustness of non-adiabatic holonomic gates 
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The robustness to different sources of error of the scheme for non-adiabatic holonomic gates 
proposed in [New J. Phys. 14, 103035 (2012)] is investigated. Open system effects as well as 
errors in the driving fields are considered. It is found that the gates can be made error resilient 
by using sufficiently short pulses. The principal limit of how short the pulses can be made is given 
by the breakdown of the quasi-monochromatic approximation. A comparison with the resilience of 
adiabatic gates is carried out. 
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I. INTRODUCTION 

A central challenge for the realization of practically 
useful quantum computation is to find systems where the 
implementation of qubits and gate operations is resilient 
to perturbations due to instabilities in the setup and to 
open system effects. Several approaches to realize this re- 
silience of quantum gates have been proposed, including 
the use of decoherence-free subspaces [T] , noisless subsys- 
tems ^ , and topological quantum computation [3] . One 
of these approaches is holonomic quantum computation 
(HQC), using adiabatic evolution [4]. This is a general 
procedure to build universal sets of gates that are ro- 
bust to certain kinds of parametric noise [5], by using 
non-Abelian adiabatic geometric phases. 

One of the most studied implementations of adiabatic 
HQC is that of Duan et at. [6], which utilizes an array 
of trapped ions that can be manipulated by laser fields. 
Schemes to implement adiabatic holonomic gates have 
also been proposed in superconducting nanocircuits us- 
ing Josephson junctions [7], and in semiconductor quan- 
tum dots 0. The gates in Refs. have turned out 
to be difficult to realize experimentally. One reason for 
this is the long run-time required for the desired para- 
metric control associated with adiabatic evolution. The 
run-time must be long enough to minimize non-adiabatic 
corrections, but sufficiently short so that the error rate 
due to open system effects and instabilities of the setup 
is small. 

A scheme for HQC based on non-adiabatic non- 
Abelian geometric phases has been proposed in Ref. [9]. 
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This scheme allows for universal quantum computation 
using gates that can be implemented rapidly compared to 
the adiabatic schemes, and therefore avoids the problems 
associated with a long run-time. Non-adiabatic HQC has 
been combined with decoherence free subspaces [10] and 
applied to coupled quantum dots as well as molecular 
magnets [ll|. In this paper, we study some aspects of 
the robustness of one-qubit gates of the proposed non- 
adiabatic scheme to different kinds of errors and compare 
it to the robustness of the corresponding adiabatic gates 
based on the Duan et al scheme. 

To realize quantum computation, it is necessary to 
bring the error rate per gate operation down below some 
threshold where error correction can be used to make the 
computation reliable. This threshold error rate has been 
estimated in different settings and by different authors 
to be somewhere between 10~^ and 10~^ [12-16]. There 
are many different sources of imperfections and their rel- 
ative importance is specific to which system is used for 
the implementation of the gate. We therefore limit our 
study to some general sources of error that are typically 
encountered in a variety of implementations. Specifically, 
we study the sensitivity to decay and dephasing and to 
imperfect control of the external driving fields. 

Previously, it has been shown that adiabatic holonomic 
quantum gates are robust to first order against small ran- 
dom perturbations of the path in parameter space [5[. 
Further analysis of robustness to parametric noise has 
been carried out in Refs. [i3-I2Q|, and robustness to envi- 
ronmental interaction has been studied in Refs. [2ll - [28j . 
Moreover, the adiabatic holonomic gates are insensitive 
to the rate at which the evolution is driven as long as the 
adiabatic approximation is valid [5]. 

The outhne of the paper is as follows. In Sec. [Ill 
we briefly review non-adiabatic and adiabatic HQC pro- 
posed in Refs. [9] and [6], respectively. In Sec. ITIH 
we study the resilience of one-qubit gates to decay of 
the excited state, dephasing, detuning, and incorrect pa- 
rameters of the driving fields. The paper ends with the 
conclusions. 
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FIG. 1: (Color online) Setup for the non-adiabatic single-qubit holonomic gate in a A configuration (left) and for the adiabatic 
single-qubit holonomic gate in a tripod configuration (right), uje is the energy spacing of the \j) and |e) level, and the driving 
fields are (jOj{t) = Vl(t)ujje^^^^ for j = 0, 1 in the non-adiabatic case and ujj{t) — Q.ujj(t/T)e^^^^ for j = 0, l,a in the adiabatic 
case. The requirement for ideal gate implementation is fj = ujje in the non-adiabatic case and z^o — cjoe = z^i — ooie = i^a — ojae 
in the adiabatic case. In the non-adiabatic setup the qubit is encoded in the |0) and |1) states. In the adiabtic setup the qubit 
is encoded in the dark subspace of the Hamiltonian. Initially the driving fields must therefore be chosen such that the dark 
subspace coincides with the computational subspace spanned by |0) and 



II. HOLONOMIC GATES 
A. The non-adiabatic gate 

The one-qubit non-adiabatic gate can be implemented 
in a A-system in the following way. A pair of zero- 
detuned pulses with the same real- valued pulse envelope 
Vt{t) couples selectively two ground state levels |0) and |1) 
to an excited state |e). The corresponding Hamiltonian 
can be expressed in the rotating wave approximation as 

= l^(i)(wo|e)(0|+wi|e)(l|+h.c.). (1) 
Here, the complex- valued driving parameters uoq and uoi 



satisfy |a;o| 



\UJi 



|2 



1, and describe the relative strength 



and relative phase of the |0) ^ |e) and |1) ^ |e) tran- 
sitions. The Hamiltonian is turned on and off at t = 
and t = r, respectively, controlled by Vt{t). The pulse en- 
velopes are described as monochromatic. It is therefore 
assumed that rt{t) is slowly varying on the time scales 
where Vj is the frequency of the driving field addressing 
the I j) ^ |e) transition, so that the quasi- monochromatic 
approximation is valid. This is equivalent to the require- 
ment that <^ 1, where /S.Vj is the spectral width of 
the pulse poj. For the rotating wave approximation to 
be valid, it is also necessary that Vt{t) <C fj [30]. Further- 
more, we take |0) and |1) to define the one-qubit state 
space. The A-configuration is illustrated in Fig. [H 

The pulse-pairs are chosen such that ujq and uji are 
time independent over the duration of each pulse pair. 
With this choice of parameters, the dark state \d) = 
—uji\0) +a;o|l) decouples from the dynamics and the evo- 
lution is reduced to a simple Rabi oscillation between the 
bright state \b) = cJq|0) +cjJ|1) and the excited state |e) 
[31|]. The Rabi frequency is ^}{t) and the subspace M{t) 
spanned by |V^^-(t)) = /o ^^"'^(*')^*' | j) = U{t,0)\j), 



j = 0, 1, undergoes a cyclic evolution if the pulse envelope 
satisfies fl{f)df = tt. Under the above conditions, the 
final time evolution operator /7(r, 0), projected onto the 
qubit space spanned by {|0), |1)}, defines the traceless 
Hermitian holonomic one-qubit gate 



U{Cn) 



n • (T, 



(2) 



where n = (sin cos ^, sin sin 0, cos 0) is a unit vec- 
tor, Cn is the evolution corresponding to uq/uji = 
— tan(^/2)e*'^, and a = {(ix^cry^az) is a vector of the 
standard Pauli operators acting on |0), Thus, any 
traceless Hermitian SU(2) operation can be implemented 
using a pulse pair. Two pairs of pulses corresponding to 
the unit vectors n and m applied sequentially results in 



U{C) = U{Cm)U{Cr.) 

= m • n + zcr • (m x n). 



(3) 



which is an arbitrary SU(2) operation. Thus, U{C) is 
a universal one-qubit gate. The evolution is purely geo- 
metric since (V^^-(t)|i^("^)(t)|V^/e(t)) = (j|i^("^)(t)|A:) = 0, 
j, /c = 0, 1, for t G [0, r]. C can be interpreted as the 
path of M{t) in the space of all 2-dimensional subspaces 
of the 3-dimensional Hilbert space, i.e., the Grassmann 
manifold G(3;2). 



B. The adiabatic gate 

The adiabatic scheme proposed by Duan et al. is 
implemented utilizing a tripod-type system, where three 
ground states |0), and \a) are coupled to an excited 
state |e) by the driving fields. Thus, the adiabatic im- 
plementation requires coherent control over an additional 
level compared to the non-adiabatic implementation us- 
ing a A-system. The tripod configuration contains a 2- 
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dimensional dark subspace, dependent on the field cou- 
plings. The system is prepared in a dark state belonging 
to the computational subspace spanned by |0) and |1), 
and the field couplings are varied independently such that 
the system remains approximately in an instantaneous 
dark state in the limit of large run-time T. More pre- 
cisely, to remain in the adiabatic regime, the evolution 
in parameter space must be slow compared to the dy- 
namical time scale of the system given by the energy dif- 
ference between the dark and bright energy eigenstates. 
The Hamiltonian for this system in the rotating wave 
approximation is 

H^-\t) = Q[cJo{t/T)\e){0\^u;,{t/T)\e){l\ 

+a;,(t/T)|e)(a| + h.c.], (4) 

where we have assumed zero-detuned field couplings. For 
ideal gate implementation it is however only necessary 
that the field couplings are equally detuned. The tripod 
configuration is illustrated in Fig. [H The field parame- 
ters are varied through a closed loop in parameter space. 
In this way, a universal set of one-qubit gates consisting 
of Ui{Ti) = e^i^ili)^!! and /72(r2) = e'^^^y , Fi and F2 
being real numbers, can be realized. 

The Ui gate can be implemented by choosing 6^0 = 0, 
ooi = — sin(i?/2)e*'^, and Ua = cos('i^/2) and taking 
and (f around a loop ci in parameter space enclosing 
the solid angle Fi = J^^sini^di^dcp. The U2 gate can 
be implemented by choosing ujq = sin'i^cos(/9, uji = 
sin sin and uja = cos'd and taking i9 and cp around 
a loop C2 in parameter space enclosing the solid angle 
F2 = J^^ sinMMip. 

III. ANALYSIS OF GATE ROBUSTNESS 

There are two qualitatively different sources of error for 
the holonomic gates. One is open system effects caused 
by the interaction with the environment. The other is 
imperfect control of the parameters of the Hamiltonian 
due to unavoidable instabilities of the setup. These pa- 
rameters are associated with the external driving fields 
used to control the system and are described as classical 
fields. 

Among open system effects, decay of the excited state 
and dephasing are important. These will be considered 
in Sec. ()IIIBp . The non-adiabatic scheme populates the 
excited state during the action of each pulse pair and it 
is therefore crucial to study the resilience to decay as a 
function of the operation parameters. In the adiabatic 
implementation, the excited state is populated only in- 
finitesimally and therefore the adiabatic gates are robust 
against decay in the adiabatic limit. While decay is typi- 
cally the dominant open system effect in trapped atomic 
systems [sS], dephasing plays a central role in supercon- 
ducting circuits [33, 34]. 

Important sources of errors due to parametric instabil- 
ity are detuning and imperfect control of the parameters 



of the driving field. These will be considered in Sec. 
pil Cp . Errors in detuning and pulse area in the non- 
adiabatic Abelian case for a two-level system have been 
analyzed in Ref. [35^. 

The gate performance under open-system effects and 
parameter errors is quantified in terms of the gate fi- 
delity (xlU'^iQgontUiQlx), where U{C) is the desired 
gate operation and ^out is the output state computed 
from the dynamics for the non-adiabatic and adiabatic 
gates with error sources. For the open system effects and 
detuning, the gate fidelities are computed numerically 
for 4000 input states |x), uniformly distributed over the 
Bloch sphere with respect to the Haar measure. The evo- 
lution of each input state is numerically integrated from 
the dynamical equations using the adaptive time step 
fourth order Runge-Kutta-Fehlberg (RKF45) method. In 
the numerical simulations the robustness to error sources 
is investigated using two test gates, the ^ phase-shift 
gate \j) j = 0,1, and the Hadamard gate 

\j) -i=[(_i)i|j) + |j0i)]^ j = 0,1. Non-adiabatic and 

adiabatic implementations of these gates are described in 

Sec. (iiira . 



A. Test gates 

1. Non-adiabatic test gates 

In the non-adiabatic scheme, the ^ phase-shift gate 
can be implemented by two pulse pairs, with the choice 
n = (cos ^, sin ^, 0) and m = (cos sin 0). The 
Hadamard gate, can be implemented by a single pulse 
pair with n = -^{^^ 0, 1). 

The TT-pulses used in the numerical simulation are hy- 
perboHc secant pulses with maximal amplitude /3. Ex- 
plicitly, for the ^ phase-shift gate we choose the two 
pulse pairs as ft{t) {ujoj'^i) = /3sech(/3t)(— 1, 1)/a/2 and 
n{t-ts)H,u[) = /3sech[/3(t - t,)](-l,e-^"/4)/V2, 
where ts is their temporal separation. For a given /3, 
the temporal separation ts must be chosen large enough 
so that the pulse overlap is negligible. This is a nec- 
essary condition to avoid any spurious dynamical con- 
tributions to the gate. The Hadamard gate is real- 
ized by a single pulse pair with shape Q{t) {ujo^uji) = 

PsechiPt) (1, {V2 - 1)) /^2{2-V2). 

The pulses are truncated where the amplitude 
is /3/1000, which gives the pulses a length r = 
2-^arcsech(Y^). Therefore, the pulse duration in the 
non-adiabatic setting decreases with increasing /3 as a re- 
sult of the pulse area being set to the fixed value tt. This 
means that the total time tr between preparation and 
read-out can be decreased as well. Thus, by increasing 
/3 we effectively speed up the gate. The implementation 
of two pulse pairs and the relevant operation times r, tg, 
and tr are illustrated in Fig. [2j Note that r < ts in order 
to avoid overlap of the pulses. 
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FIG. 2: (Color online) Schematic illustration of the imple- 
mentation of a non-adiabatic holonomic gate. The time tr 
between preparation (P) and read-out (R) of the qubit, the 
separation ts of the pulse-pairs and the pulse length r are 
indicated. 




(i?, ^) = (0, 0) ^ (f , 0) ^ (f, ^) ^ (0, tt) ^ (0, 0) 




{i}, ^) = (0, 0) ^ (f , 0) ^ ip) = (0, 0) (f , 0) ^ 
(I. -!) - (0, -!) - (0, 0) (f, 2^) ^ (0, 2^) ^ (0, 0) 

FIG. 3: (Color online) The curve in parameter space used to 
implement the adiabatic holonomic |- phase-shift gate (top) 
and the two curves in parameter space used to implement the 
Hadamard gate (bottom). The top and right bottom curves 
correspond to the field couplings cjq = 0, cji = — sin (1^/2) e*'^ 
and uoa = cos('i^/2), while the left bottom curve corresponds 
to (jJo = sin-i^ cos cji = sin 1^ sin and uja = cosi^. 



2. Adiabatic test gates 

Next, we consider the adiabatic scheme. The ^ phase- 
shift gate can be implemented by one loop in parame- 
ter space, while the Hadamard gate requires two loops. 
The ideal adiabatic ^ phase-shift gate is generated in 
the 00 Hmit by varying the field couplings cjq = 0, 

uji = — sin(i?/2)e*'^ and uJa = cos('i^/2) along the loop 
(^,(^) = (0,0) ^ (|,0) ^ (f,^) ^ (0,^) ^ (0,0) at 
constant speed. Here, |0) is decoupled from the excited 
state. 

The ideal Hadamard gate is generated in the QT 00 



limit by combining two loops. In the first loop, the field 
couplings are ujq = sin'i^cos(/9, uji = sini^smcp and uJa = 
cos i9 and these are varied along the loop {'d^ (p) = (0, 0) 
(f,0) ^ (|,-f) (0,-f) ^ (0,0) at constant speed. 
In the second loop, the field couplings are ujq = 0, cji = 
— sin(i^/2)e*'^ and uJa = cos(i?/2) and these are varied 
along the loop {i9,(p) = (0,0) ^ (f,0) ^ (f,27r) ^ 
(0,27r) (0,0). The curves in parameter space used to 
implement the adiabatic gates are illustrated in Fig. [3l 

Since the run-time T and the strength of the field cou- 
pling Q are always finite for adiabatic gates, there will be 
non-adiabatic corrections that reduce the fidelity. These 
corrections are due to failure of the state to remain in 
the dark subspace during evolution. In Fig. HI we show 
the fidelity for the adiabatic test gates, for finite Q and 
T. The fidelities are plotted as functions of the dimen- 
sionless quantity QT. 

The oscillatory behavior of the fidelities, as a function 
of l^T, is due to non-adiabatic effects. These effects can 
be understood as a combination of modified holonomies, 
caused by changes in the path of the computational sub- 
space, and non-zero dynamical phases. The revivals seen 
in the fidelities have been pointed out previously in Ref. 
[25| . Some of the revivals reach unit fidelity and these 
can therefore be used to implement the desired gates for 
finite QT. However in the non-adiabatic regime the gate 
is not holonomic since there are both dynamical and ge- 
ometrical contributions to the gate operation. For the ^ 
phase-shift gate the maximal fidelity (red line) is unity 
for all since |0) is decoupled from the dynamics. 

B. Open system effects 

1. Decay 

To investigate the sensitivity of the non-adiabatic 
scheme to decay we assume that the time scale of the 
applied pulses is large compared to the time scale of 
the dynamics underlying the decay process, so that the 
Markovian approximation is valid. Given this, we fur- 
ther assume that the excited state decays to an auxil- 
iary ground state level \g) js^ with a time independent 
rate 7. We compare the sensitivity of the non-adiabatic 
implementations of the test gates with that of the cor- 
responding adiabatic implementations of the gates. The 
decay is modelled by the Lindblad equation 

Qt = -i[H{t), Qt] + 2LgtL^ - L^Lgt - QtL^L, (5) 

where Qt is the density operator, L = y^|^)(e|, and H(t) 
is either i^^"^)(t) or H^^'^t). 

In the non-adiabatic case, the constraint of having tt- 
pulses implies that the only experimentally controllable 
parameters of principal importance are the maximal cou- 
pling strength /3 and the total time between preparation 
and read-out t^. The dynamics of the gate can be de- 
scribed in terms of two dimensionless parameters. These 
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FIG. 4: (Color online) Influence of non-adiabatic corrections on the adiabatic holonomic ^ phase-shift gate (left) and Hadamard 
gate (right). The effect is quantified from top to bottom in terms of maximum (red), average (black), and minimum (blue) 
fidelities. We plot the fidelities as functions of the dimensionless quantity QT, where Q is the time independent global strength 
of field couplings and T is the run-time of the gate. 



can be chosen as /3/7 and jU. When the pulse is a per- 
fect TT-pulse the excited state is negligibly populated after 
the end of the pulse. Therefore, the relevant time scale 
for decay is the width of the pulse r, which is propor- 
tional to ^. Since jr (x ^ there is only one dynamically 
relevant dimensionless parameter. 

The relevant operation parameters in the adiabatic 
case are the coupling strength and run-time T. The 
dynamics can be characterized by two dimensionless pa- 
rameters, which can be chosen as QT and 7T. In the 
simulations we consider both a fixed coupling strength 
and vary T, and a fixed run- time Tq and vary Q. 

In Fig. [5l we show the fidelities of the test gates, 
computed from Eq. (|5|). The fidelities are plotted as 
functions of the dimensionless quantity /3/7 in the non- 
adiabatic case, as well as ftoT and ^ITq in the adia- 
batic case for fix coupling strength and fixed run- time, 
respectively. For the adiabatic case, we have chosen 
flo/j = 12.5, in the case of fixed coupling strength, and 
in the case of fixed run-time we have chosen 7T0 = 8 for 
the ^ phase-shift gate and 7T0 = 32 for the Hadamard 
gate. These choices are made to make the effect of decay 
non- negligible. Furthermore, in the non-adiabatic case 
we have chosen 7^^ = 8, which guarantees that there is 
no pulse overlap for the f3/j range shown. 

The fidelities of the non-adiabatic gates tend monoton- 
ically to unity in the large P/j limit for both test gates 
(left panels). This demonstrates that the non-adiabatic 
holonomic test gates can be made robust to decay of the 
excited state by employing pulses that are sufficiently 
short compared to the decay time ^ . For the Hadamard 
gate the maximal fidelity (red line) is identity for all /3/7. 
This is because it is implemented by a single pulse and 
the dark state of the Hamiltonian corresponding to this 
pulse is left unchanged by the dynamics and hence unaf- 
fected by the decay. For the ^ phase-shift gate no state 
is left unchanged by the gate operation, which explains 
why maximal fidelity is slightly decreased for small f3/j. 

The stability of the adiabatic gates to decay of the ex- 
cited state in the adiabatic limit is confirmed as the fideli- 



ties of the adiabatic gates tend to unity in the large QT 
limit. However, the fidelity of the revivals will no longer 
reach unity when there is decay of the excited state. In 
the adiabatic implementation of the ^ phase-shift gate 
the maximal fidelity (red line) is identity since |0) is de- 
coupled from the dynamics and therefore unaffected by 
the decay. For the Hadamard gate no state is decoupled 
from the dynamics and therefore the maximum fidelity 
is low for small QT. 



2. Dephasing 

Dephasing is hard to eliminate in some implementa- 
tions of holonomic gates, for example in superconducting 
Josephson junctions [3^. Therefore, we study the ro- 
bustness of the non-adiabatic and adiabatic schemes to 
dephasing. More precisely, we consider dephasing in the 
|0), |e) and |e) bases. The effect of dephasing is mod- 
eled by the Lindblad equation 



Qt 



( 



i=0,l 



2LkeQtL 



LlLi 



■Qt 



^,4L,e),(6) 



where Qt is the density operator, Lke = V^(|e)(e| — \k){k\) 
are the Lindblad operators, and H(t) is either H'^^^^t) 
or H^^\t). We use the assumption that the time scale of 
the process underlying the dephasing is short compared 
to the time scale of the pulses (Markovian approxima- 
tion). We again compare the resulting fidelities of the 
non-adiabatic implementations of the two test gates with 
those of their corresponding adiabatic implementations. 

In the non-adiabatic case, the operation parameters 
are P and the total time tr between preparation and 
read-out. The two dimensionless parameters describing 
the dynamics are | and et^. Since the qubit space is 
not a decoherence free subspace of the dephasing, the 
full time between preparation and read-out is relevant 
for the fidelity of the gate operation. The best fidelity 
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FIG. 5: (Color online) Influence of decay with rate 7 of the excited state |e) on the non-adiabatic and adiabatic holonomic ^ 
phase-shift gate (upper) and Hadamard gate (lower). The effect is quantified from top to bottom in terms of maximum (red), 
average (black), and minimum (blue) fidelities. The three panels show from left to right, the non-adiabatic gate with decay 
(left) and the adiabatic gate with decay, for a fixed coupling strength (middle) and fixed run-time (right). Choosing hyperbolic 
secant 7r-pulses with amplitude /3, the non-adiabatic fidelities are plotted as functions of the dimensionless quantity P/j. We 
plot the adiabatic fidelities as functions of the dimensionless quantities QqT and QTo, where Q is the time independent global 
strength of field couplings, T is the run-time of the gate, and Qo and To are particular fixed values of these quantities. For 
the adiabatic gates we have chosen ^0/7 = 12.5 for the case with fixed coupling strength and for the case of fixed run-time 
we have chosen 7T0 = 8 for the ^ phase-shift gate and 7T0 = 32 for the Hadamard gate. In the non-adiabatic case, we have 
chosen jts = 8, where ts is the temporal separation of the two pulses, ts is chosen sufficiently large to avoid pulse overlap for 
the /3/7 range shown in the left panels. 



is thus achieved when the total time is reduced to only 
the time required to implement the gates. To study the 
effect of dephasing on gate operation we have therefore 
assumed that gates are implemented immediately follow- 
ing preparation and that read-out is made immediately 
afterwards. This assumption reduces the relevant time 
to r cx and the only relevant dimensionless parameter 

is therefore -. 

e 

The relevant operation parameters in the adiabatic 
case are the coupHng strength Q and run- time T, and the 
dimensionless parameters describing the dynamics can be 
chosen as QT and eT. In the simulations we consider both 
a fixed coupling strength and vary the run- time, and 
a fixed run-time Tq and vary the coupling strength. 

In Fig. [6l we show the fidelities of the test gates, 
computed using Eq. (j6j). The fidelities are plotted as 
functions of the dimensionless quantity /3/e in the non- 
adiabatic case, as well as rtoT and QTq in the adiabatic 
case for fix coupling strength and fixed run- time, respec- 
tively. For the adiabatic case we have chosen fto/e = 
78.125, in the case of fixed coupling strength, and in the 
case of fixed run- time we have chosen cTq = 0.16 for the 
^ phase-shift gate and cTq = 0.32 for the Hadamard gate. 
These choices are made to make the effect of dephasing 
non- negligible. Furthermore, in the non-adiabatic case 
we have chosen ets = 0.128, which guarantees that there 
is no pulse overlap for the /3/e range shown. 



The fidelities of the non-adiabatic gates tend mono- 
tonically to unity in the large P/e limit for both test 
gates (left panels). Thus, the non-adiabatic version of 
the holonomic ^ phase-shift gate can be made resilient 
to the dephasing channel by employing sufficiently short 
pulses and reducing idle time before and after the pulses. 

The adiabatic gates are not stable to dephasing in the 
limit T ^ 00 and we see a gradual decline in the fidelities 
as run- time increases. If the run- time is fixed and Q 
is increased the fidelities stabilize at some value below 
unity that is a function of the parameter cTq. To have 
high fidelity in the adiabatic case it is necessary that Q 
can be made large compared to e so that the adiabatic 
approximation is valid for a run-time at which the effect 
of dephasing is still negligible. Alternatively, Cl should be 
large enough so that the effect of dephasing is negligible 
at some T corresponding to a revival of the fidelity. 



C. Parametric control 

1. Detuning 

We assume that the \j) ^ \e) transition is driven by a 
laser pulse with frequency uj. The associated detuning 
is Aj = 27TUj — ujje^ where ujje is the corresponding en- 
ergy spacing. Ideally, the detunings must vanish in the 



7 




FIG. 6: (Color online) Influence of dephasing in the |0), |e) and |1), |e) bases, on the non-adiabatic and adiabatic holonomic ^ 
phase-shift gate (upper) and Hadamard gate (lower). The effect is quantified from top to bottom in terms of maximum (red), 
average (black), and minimum (blue) fidelities. The three panels show from left to right, the non-adiabatic gate with dephasing 
(left) and the adiabatic gate with dephasing, for a fixed coupling strength (middle) and fixed run-time (right). Choosing 
hyperbolic secant 7r-pulses with amplitude /3, the non-adiabatic fidelities are plotted as functions of the dimensionless quantity 
P/e. We plot the adiabatic fidelities as functions of the dimensionless quantities QqT and QTo, where Q is the time independent 
global strength of field couplings, T is the run-time of the gate, and Qo and To are particular fixed values of these quantities. 
For the adiabatic gates, we have chosen Qo/e = 78.125 for the case with fixed coupling strength and for the case of fixed 
run-time we have chosen eTo = 0.16 for the f phase-shift gate and eTo = 0.32 for the Hadamard gate. In the non-adiabatic 
case, we have chosen ets = 0.128, where ts is the temporal separation of the two pulses, ts is chosen sufficiently large to avoid 
pulse overlap for the /3/e range shown in the left panels. 



non-adiabatic case while they must be all equal in the 
adiabatic setting. Here, we examine the effect of devi- 
ations from these ideal values on the gate fidelity. To 
simplify the analysis, we limit the study to time inde- 
pendent detunings. 

Non-zero detunings give rise to additional diagonal 
terms in the Hamiltonian. In the non-adiabatic case, we 
have 

H^r\t) = Ao|0)(0|+Ai|l)(l|+J?("-)(t), (7) 

where H^^^\t) is the ideal Hamiltonian in Eq. ([1]). Sim- 
ilarly, in the adiabatic case, we have 

H^^\t) = Ao|0){0| + Ai|l){l| + A„|a){a|+ff(-)(t), 

(8) 

where now H'^^^t) is the ideal Hamiltonian in Eq. 
If two detunings are equal but different from the third, 
there will only be one dark state, and if all three detun- 
ings are different from each other there will be no dark 
state at all. Thus, deviations from the Aq = Ai = A^ 
constraint destroy the dark state structure and the asso- 
ciated holonomy. 

Note that we express H^^\t) and H^\t) in frames 
that are "co-rotating" with their respective detuned driv- 
ing fields. Therefore, to compare the output of the 



gate operation ^out(^) generated by the detuned Hamil- 
tonians, with the output of the desired gate operation 
U (C), we must transform ^out(^) to the frame co-rotating 
with the ideal driving fields. This is done through 
the transformation ^out(^) e*'^^^out(^)e~*'^S where 
S = Ao|0)(0| + Ai|l)(l| in the non-adiabatic case and 
S = Ao|0)(0| + Ai|l)(l| -I- Aa|a)(a| in the adiabatic case. 
Note also that these frame rotations will remove the ef- 
fect of the diagonal terms in Eqs. ([7]) and ([8]) for the part 
of the evolution where the driving fields vanish. 

We study two principal cases. First, all detunings are 
set equal, and secondly, some of the detunings are as- 
sumed to be different. These cases are naturally captured 
by the mean and relative detunings. In the non-adiabatic 
setting, these read A = and 6 = , respec- 

tively. In the adiabatic case, we have the mean detuning 
Aoia = ^"^^^^^"^ two independent relative detun- 
ings 5oi = and 5oa = 2^"" » In some implemen- 
tations it may be easier to control the relative detuning 
of the driving fields, than the mean detuning. 

First, we investigate how gate operation is affected in 
the non-adiabatic case when a constant mean detuning is 
introduced and the relative detuning is zero. For compar- 
ison we also include the adiabatic implementations with 
the same mean detuning and all relative detunings zero. 

The dynamics of the non-adiabatic gate with A 7^ 
and S = depends on two dimensionless parameters that 
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FIG. 7: (Color online) Influence of a constant mean detuning A and zero relative detuning on the non-adiabatic and adiabatic 
holonomic ^ phase-shift gate (upper) and Hadamard gate (lower). The effect is quantified from top to bottom in terms of 
maximum (red), average (black), and minimum (blue) fidelities. The three panels show from left to right, the non-adiabatic 
gate with mean detuning (left) and the adiabatic gate with mean detuning, for a fixed coupling strength (middle) and fixed 
run-time (right). Choosing hyperbolic secant vr-pulses with amplitude /3, the non-adiabatic fidelities are plotted as functions 
of the dimensionless quantity 13/ A. We plot the adiabatic fidelities as functions of the dimensionless quantities QqT and QTo, 
where Q is the time independent global strength of field couplings, T is the run-time of the gate, and Qo and To are particular 
fixed values of these quantities. For the adiabatic gates we have chosen Qo/Aoia = 6.25 for the case with fixed coupling strength 
and for the case of fixed run-time we have chosen AoiaTo = 16 for the |- phase-shift gate and AoiaTb = 64 for the Hadamard 
gate. In the non-adiabatic case we have chosen Ats = 80, where ts is the temporal separation of the two pulses, ts is chosen 
sufficiently large to avoid pulse overlap for the 13/ A range shown in the left panels. 



can be chosen as /3/A and Atr. However since the entire 
dynamics is generated by the driving fields the relevant 
time is the pulse length r oc ^. Thus, /3/A is the only 
relevant dimensionless parameter. 

The relevant operation parameters in the adiabatic 
case are the coupling strength Q and run-time T, and 
the dimensionless parameters describing the dynamics 
are QT and AT. In the simulations we consider both 
a fixed coupling strength Qq and vary the run-time and 
a fixed run-time To and vary the coupling strength. 

In Fig. [71 we show the fidelities of the gates com- 
puted using Eqs. ([7j) and ([8|) in the non-adiabatic and 
adiabatic case, respectively. The fidelities are plotted as 
functions of the dimensionless quantity /3/A in the non- 
adiabatic case, as well as ftoT and QTq in the adiabatic 
case for fix coupling strength and fixed run- time, respec- 
tively. We choose l^o/Aoia = 6.25 in the case of fixed 
coupling strength. In the case of fixed run-time, we take 
^oia^o = 16 for the ^ phase-shift gate and Aoia^o = 64 
for the Hadamard gate. Furthermore Ats = 80, which 
guarantees that there is no pulse overlap for the P/A 
range shown in the figure. 

With a constant mean detuning and zero relative de- 
tuning, the fidelities of the non-adiabatic gates tend to 
unity in the large P/A limit (left panels). This demon- 
strates that the non-adiabatic versions of the holonomic 
test gates can be made robust to constant mean detun- 



ing by employing sufficiently short pulses. The adia- 
batic scheme is stable to constant mean detuning in the 
QT oo limit. 

Next, we consider nonzero relative detuning. The adi- 
abatic implementation of the ^ phase-shift gate does not 
involve any driving field coupled to the |0) ^ |e) transi- 
tion. For this reason we consider the case where the only 
nonzero detuning is (5 = Ai of the driving field associated 
with the |1) ^ |e) transition. 

The relevant dimensionless dynamical parameter for 
the non-adiabatic gates is /3/6. The relevant opera- 
tion parameters in the adiabatic case are the coupling 
strength Q and run- time T, and the dimensionless pa- 
rameters describing the dynamics are QT and 5T. In the 
simulations we consider both a fixed coupling strength 
Qq and vary T and a fixed run-time Tq and vary Q. 

In Fig. [HI we show the fidelities of the test gates com- 
puted using Eq. ([7j) and ([8]) in the non-adiabatic and 
adiabatic case, respectively. The fidelities are plotted as 
functions of the dimensionless quantity P/S in the non- 
adiabatic case, as well as ftoT and QTq in the adiabatic 
case for fix coupling strength and fixed run- time, respec- 
tively. We have chosen Qq/S = 6.25 in the case of fixed 
coupHng strength, while in the case of fixed run-time, we 
take STq = 2 for the ^ phase-shift gate and STq = 64 for 
the Hadamard gate. Furthermore, Stg = 16, which guar- 
antees that there is no pulse overlap for the /3/d range 
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FIG. 8: (Color online) Influence of a constant relative detuning 6 on the non-adiabatic and adiabatic holonomic |- phase-shift 
gate (upper) and Hadamard gate (lower). The effect is quantified from top to bottom in terms of maximum (red), average 
(black), and minimum (blue) fidelities. The three panels show, from left to right, the non-adiabatic gate with relative detuning 
and the adiabatic gate with relative detuning, for a fixed coupling strength (middle) and fixed run-time (right). Choosing 
hyperbolic secant vr-pulses with amplitude /3, the non-adiabatic fidelities are plotted as functions of the dimensionless quantity 
13/6. We plot the adiabatic fidelities as functions of the dimensionless quantities QqT and QTo, where Q is the time independent 
global strength of field couplings, T is the run-time of the gate, and Qo and To are particular fixed values of these quantities. 
For the adiabatic gates we have chosen Qo/6 = 6.25 for the case with fixed coupling strength and for the case of fixed run-time 
we have chosen 6To = 2 for the ^ phase-shift gate and 6To = 64 for the Hadamard gate. In the non-adiabatic case we have 
chosen Sts = 16, where ts is the temporal separation of the two pulses, ts is chosen sufficiently large to avoid pulse overlap for 
the /3/6 range shown in the left panels. 



shown in the figure. 



With constant relative detuning, the fidelities of the 
non-adiabatic gates tend to unity in the large P/S limit 
(left panels). The behavior is similar to the case with 
nonzero mean detuning. The adiabatic gates, on the 
other hand, are unstable to relative detuning and do not 
converge to any value of the fidelity when run-time T is 
increased while the coupling strength is fixed. Instead, 
the mean fidelity as a function of T is an oscillating func- 
tion. If the run-time is fixed and ft is increased the fi- 
delities stabilize at some value that is a function of the 
parameter STq and typically not unity. 



The non-adiabatic gates can thus be made resilient to 
relative detuning by using short enough pulses. One way 
to have high fidelity in the adiabatic case, is to choose 
T corresponding to a maximum in the oscillating mean 
fidelity. This however requires precise knowledge of the 
relative detuning. Without such knowledge it must be 
possible to choose ft/ 6 sufficiently large so that the first 
decline in fidelity due to relative detuning becomes sig- 
nificant only at run-times larger than the T required for 
the adiabatic approximation to be valid. 



2. Driving field parameters 

Control of the pulse envelope and the relative ampli- 
tudes and phases of the two driving fields is of crucial 
importance to gate operation. To study the effect of er- 
rors in the driving fields, we make the assumption that 
the deviations from the ideal case are such that the rel- 
ative strength and relative phase of the |0) ^ |e) and 
|1) ^ |e) transitions are time independent during the 
implementation of a pulse-pair. Given this assumption 
the Hamiltonian can be written 

H{t) = l^(t)(cI}o|e)(0|+Sii|e)(l|+h.c.) 

= m^o. (9) 

where (l{t) is the non-ideal pulse envelope and ujj are the 
non-ideal relative strength and phase of the transitions 
satisfying |cDoP + l^iP = 1- Since the implementation of 
the non-adiabatic gates depends only on the area of the 
two pulses in each pulse pair and not on the exact shape, 
gate operation is robust to area preserving deviations in 
the shape. Only the deviation of a = (l{t)dt from tt is 
relevant. 

If we introduce the notation H = Q{t)Ho for the ideal 
Hamiltonian we can express the error due to the devia- 
tions in terms of the fidelity as 

= |(V^|e^"^Oe-^^^°(^)|^)p, (10) 
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where is an arbitrary normalized input qubit state. 
Using that i^^^ = and HI""-^ = Hq for n = 1^2, . . . , 
we can see that 6^^^° = 1 - 2H^. Similarly H^"^ = 
Hi and ^0^^^^ = implies that e-^^^^ = !-(!- 
cosa)i^Q — z(sina)i^o. Using this and the fact that the 
expectation value with respect to |?/^) of any product of 
an odd number of Hq and Hq vanishes, we can express 
the fidelity as 



\1-2m 



-(1 



COS 



a) (. 



(11) 



Assuming that a = ir -\- Sa and Hq = Hq -\- SHq^ where 
6a and SHq = ^a;o|e)(0| + (5a;i|e)(l| + h.c. are small, we 
can expand the fidelity to second order in the deviations. 
This gives 



1- Sa^{i>\H^\i. 



)-4{^\5Hi\^) 
5HlHl\^) 



(12) 



Averaging over all input states gives the average fidelity 
1 



l-±Sa^-2{\Suof + \Su,f) 



(13) 



Thus, we can see that the one-qubit non-adiabatic holo- 
nomic gate is robust to first order in the deviations in the 
pulse area as well as in the relative phase and strength 
of the transitions. The error incurred by the incorrect 
parameters can be understood as a failure of the sub- 
space M{t) to follow the correct path. If the area of the 
pulses deviates from tt the state will not return to the 
computational subspace, and the final state will have a 
nonzero amplitude in the excited state. If the pulse area 
is TT, the state will return to the computational subspace 
but the gate operation will not be the desired one unless 
— = — = e*" for some a G M amounting^ to a e*" phase 

UJQ UJi o jr 

shift of both ujQ and uJi. 

A special case is when there is a deviation in the pulse 
area but correct relative strength and phase of the cou- 
pHngs. The second order dependence of the fidelity on 
Sa in this case is consistent with the result of Ref. [35| 
for the Abelian case. 

Although the above deviations in the parameters can 
lead to a gate operation that is not the desired one, the 
evolution is nevertheless still purely geometric during the 
implementation of the gate. This is in contrast to the er- 
ror caused by detuning, where the reduced fidelity is due 
to the combined effect of modified holonomies and dy- 
namical phases, and to the case with open system effects 
where the state evolves into a mixture of states that pick 
up different combinations of holonomies and dynamical 
phases. Another difference is that the error introduced 
by incorrect parameters is independent of the run-time 
of the gate. 



IV. CONCLUSIONS 

The non-adiabatic gate can be made resilient to decay 
of the excited state and to constant mean and relative 
detunings by employing pulses that are sufficiently short 
compared to the time scales of the decay and detuning. If 
the idle time between preparation of the qubit and read- 
out can be made negligible, the gate will also be resilient 
to dephasing in the |e), |1) and |e), |0) bases, in the limit 
of short pulses. It is therefore of critical importance for 
the implementation of such gates that the pulse height /3 
can be made sufficiently large relative to the dynamical 
parameters describing these sources of error. 

There is a principal upper limit on how fast the pulses 
can be implemented given by the breakdown of the quasi- 
monochromatic approximation when the pulse changes 
rapidly compared to the oscillations of the driving field. 
If the quasi- monochromatic approximation fails there will 
be non-negligible frequency components of the driving 
field other than the desired one. These may couple to 
other transitions in the system and reduce fidelity of the 
gate. There is also a limit given by the breakdown of the 
rotating wave approximation when the ratio of the pulse 
height to either of the energy spacings ujje^ j = or 1, 
becomes too large. This causes a non negligible Bloch- 
Siegert shift j30|,[37|] of the \j) ^ \e) transition resonance 
frequency due to the effect of the counter-rotating terms 
that can no longer be neglected. Furthermore, the anal- 
ysis of the sensitivity to open-system effects is only valid 
when the Markovian approximation holds. Thus, if the 
pulse duration is decreased to a point where it becomes 
comparable to the time scale of the dynamical processes 
underlying the open system effects, the analysis using 
Lindblad's equation is no longer valid and memory ef- 
fects of the environment must be taken into account. 

Control over the pulse shape is important only to the 
extent that the pulse area must be tt. Moreover, the 
fidelity depends on small deviations in the pulse area and 
the other parameters of the driving fields only to second 
order. The gate cannot be made more robust to this 
source of error by decreasing run- time. 

In comparison, the adiabatic gates are robust to de- 
cay and mean detuning in the limit of large run- time. 
They are however not robust to dephasing and relative 
detuning in this limit. Resilience to these sources of er- 
ror in the adiabatic case requires that the field coupling 
strengths can be made large compared to the dephasing 
and relative detuning parameters, and that the run- time 
can be chosen sufficiently small to make the error due to 
relative detuning or dephasing negligible. The require- 
ments on the operation parameters for high fidelity gate 
implementation in the presence of dephasing or relative 
detuning are thus qualitatively similar to those for the 
non-adiabatic gate in the sense that the coupling strength 
must be increased and the run-time decreased. 

To fully address the issue of robustness of the non- 
adiabatic scheme for quantum computation one must also 
analyze the robustness of the two-qubit gate proposed in 
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Ref. jil. The two-qubit gate involves coupled A systems, 
that could be implemented using for example trapped 
ions coupled via collective spatial vibrations as in the 
S0rensen-M0lmer ion trap scheme [ss]. In addition to 
errors emanating from imperfections in the driving fields 
and the interaction of the A systems with the environ- 
ment, the analysis would have to involve errors originat- 
ing from the coupling mechanism as well. 
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